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Abstract 



O 

CD 

Q . We consider a one dimensional elastic string as a set of massless beads in- 



teracting through springs characterized by anisotropic elastic constants. The 



^ . string, driven by an external force, moves in a medium with quenched disor- 

^ . der. We present evidence that the consideration of longitudinal fluctuations 

leads to nonlinear behavior in the equation of motion which is kinematically 
0^ ■ generated by the motion of the string. The strength of the nonlinear effects 



depends on the anisotropy of the medium and the distance from the depinning 
transition. On the other hand the consideration of restricted solid on solid 



O ■ conditions imposed to the growth of the string leads to a nonlinear term in the 



equation of motion with a diverging coefficient at the depinning transition. 
The motion of an elastic string in disordered media has attracted considerable recent 
attention, in part due to its relevance to flux-flow in type-II superconductors [jl| and rough- 
ening of non-equilibrium interfaces 0. By means of a number of numerical and ana- 
lytical P,p!0| studies it has been observed that scaling theory can be used as an underlying 



framework to understand and characterize the dynamical properties of the elastic string. 

Consider a one- dimensional elastic string moving under the influence of an external 
driving force F, normal to the string, in a two-dimensional disordered medium of edge 
L. A discrete model for such a string consists on massless beads connected by springs. 
The string is assumed to be oriented along the x-axis and the position of the i-th bead 
is denoted by a two-dimensional displacement vector rj = {xi,yi), i = 1,...,L (see Fig. 



H]). The disorder in the medium is introduced by uniformly-distributed pinning sites with 
random strength, which we refer to as quenched disorder or "quenched noise" . The dynamics 
of such a string is the result of the interplay between the quenched disorder characteristic 
of the medium and the elastic properties of the string. 

A key quantity is the average velocity of the string as a function of the external force. At 
small forces F the string is pinned by static disorder. Just above the depinning transition 
F = Fc, i.e. when the external force overcomes the pinning effect of impurities, the velocity 



varies as 



vo{f)-f, (1) 

where 6 is the velocity exponent and / = F/F^ — 1 the reduced force. 

Neglecting thermal fluctuations and lateral fluctuations of the beads, the equation of 



motion for the string in the continuum limit is the Edwards- Wilkinson equation [|TT| with 
quenched disorder [3-|T0| 

^y^ = uV'y + vix,y) + F (2) 

The first term in the right hand side of (0) includes the elastic effects acting to make the 
string straight. The second term mimics the quenched disorder, which has zero mean and is 
uncorrelated. The string is driven in the y direction by the external force F. For large driving 
force (F ^ F^), the quenched noise becomes effectively time-dependent, ri{x,yQ + vt). It is 
believed |[TO|Jl^ that in this regime the motion of the string induces an additional nonlinear 
term in (Q), namely the Kardar-Parisi- Zhang (KPZ) term A(Vy)^ Jl^. However, since this 
nonlinear term is generated by the motion of the string, A is expected to vanish as the 
velocity goes to zero at the depinning transition, and the critical behavior at F = F^. is 
correctly described by Eq. (^. 

While d^) can be obtained (using dty = —6T-C/6y + F) from the Hamiltonian 

n= f"^ dx{v{Vyf + ^i{x,y)}, (3) 
Jo 



the KPZ nonlinear term A(V?/)^ cannot be deduced as a variation of any bounded Hamilto- 
nian. Here the quenched noise is ri{x,y) = —6yfi{x,y). 

In the Hamiltonian (|^), only transverse fluctuations (along the y direction) contribute 
to the elastic energy, forbidding longitudinal fluctuations (along the x direction). However, 
for a real elastic string, the elastic energy depends on the distance (r^ — ri_i)^ between 
two consecutives beads. Here we introduce a (l+l)-dimensional model which allows both 
for longitudinal and transverse fluctuations of the beads. In the model, the elastic energy 
depends both on Uxi^i — Xj_i)^ and I'y^yi — yi-iY, where and Uy are the elastic constants 
corresponding to displacements in the x and y direction, respectively. We focus on the 
determination of the equation of motion of the string. We flnd that, even though the 
string can form overhangs, at large enough length scales the string still have a well-defined 
orientation and profile, and can be described by a continuum theory. The main results of 
this paper are: 

(a) In the limit, e = Vy/vx ^ 1? where e is the anisotropy parameter, the large-scale 
behavior of the string is described by the nonlinear equation of motion with quenched noise 

= uV'y + KVyf + r^(x, y) + F, (4) 

where the nonlinear term A(V?/)^ in (^ is of kinematic origin. We find that A vanishes at 
the depinning transition as 

A(/) ~ /l^l - 0. (5) 

(b) If longitudinal fiuctuations are neglected, nonlinear terms of the type A(Vy)^ are for- 
bidden in the growth equation. We argue that this result applies to a number of previously- 
introduced models 0-§[. In our model, this limit corresponds to e = Uy/ux 0. 

(c) A different scenario is found when the rules of motion of the beads are constraint 



to satisfy a restricted sohd on solid (RSOS) condition \hi±i — hi\ < const |T^. When such 
condition is imposed we find that the equation of motion of the string is (^ but with a 
coefficient A which diverges at the depinning transition as 



A(/) ~ f-'^ ^ oo. (6) 
This result is valid for any value of the anisotropy parameter e, and applies to a number of 



growth model in the directed percolation universality class |[T5|-pO . 

We now take up each of these results in turn. Before beginning, we note that for a given 
model the presence of a nonlinear term A(V?/)^ can be identified using tilt-dependent velocity 



measurements pT|JT8| -pO|. Suppose we tilt the elastic string, by imposing helical boundary 
conditions yi = Vl + mL, where m is the average tilt of the string. Then, according to @, 
the average tilt-dependent velocity becomes 

v{m) = vq + Am^, (7) 

where vq is the velocity of the untilted string. If A = 0, so that the motion of the elastic string 
is described by (^, then the velocity does not depend on the average tilt of the interface. 
Tilt dependence is expected only if there is a nonlinear term in the equation of motion of the 
form A(V?/)^. This property can be used to gain information on the presence and magnitude 
of the nonlinear term A, by monitoring the velocity of the string as a function of the average 
tilt, and fitting to a parabola the obtained curve [p2| . 



In the following, we study a generalized model of an elastic string that allows for lateral 
motions of the beads and therefore overhangs. The main element of the model, not included 
in the Hamiltonian (^, is the existence of longitudinal motion of the beads. To include this 
additional degree of freedom, we use a generalized Hamiltonian 

(8) 

IJ,{xi,yi) - Fyi]. 

We simulate the discrete version of (P), concentrating on the zero-temperature dynamics 
of the string (only motions which decrease the total energy of the string are allowed). A 
standard Monte Carlo algorithm, by choosing randomly a site on the interface, induces time- 
dependent noise. Since at zero temperature the motion of the string is deterministic, we 
have employed an algorithm with parallel updating, during which even and odd sublattices 



are updated simultaneously. The quenched noise is introduced by defining at every site 
of the two-dimensional lattice uncorrelated random numbers uniformly distributed 

between —6 and 6. 

During the simulations, the chosen bead is allowed to move to one of its four nearest 
neighbors, if that motion decreases the total energy of the string given by (^. If there are 
more than one possible moves with ATi < 0, then the one with most negative ATi is chosen. 
We focus on the determination of the nonlinear term A, measuring the tilt-dependent velocity 
of the string. 

(a) Figure ^ shows the velocity of the driven elastic string as a function of the average tilt 
for different driving forces. The results correspond to the anisotropic motion characterized 
by i^a; = 0.1 and Uy = 1 {e = 10). We see that the velocity follows a parabola with the tilt, 
indicating the presence of a nonlinear term A(V?/)^ above the depinning transition (moving 
phase, F > Fc). However, the parabolas become flatter as the depinning transition is 
approached. Our calculations indicate that A ^ as F — as in (|^). These results are 
obtained for the anisotropic case z/^. < Uy, {e > 1), and further increasing the anisotropy, the 
observed behavior does not vanish. 

(b) The other limit of the model leads to known results: a flnite Uy and i/^, — > oo means 
that longitudinal fluctuations are energetically very expensive, allowing only transversal 
fluctuations. In this limit, the model reduces to the models of Refs. [^|-|^, where longitudinal 
fluctuations are not allowed. In this case the nonlinear term is exactly zero (see Fig. 
Thus as i/j: — > oo, a slow decrease of A toward zero is expected. Our simulations indicate 
that the decrease of A is much faster, dropping to immeasurably small values as we approach 
the isotropic point, z/^. = Uy. 

There are two possible scenarios compatible with these results. According to the flrst, 
A = for £ < 1 and A 7^ for e > 1. The second scenario says that A ^ as £ ^ 0. Our 
data support the flrst hypothesis, but at this point we cannot exclude the second. 

(c) Figure l^c shows the results of our simulations when the RSOS condition is applied to 
the growth of the string: for a given i-th bead, if /liii — hi>2 then we increment hi ^ hi + 1 
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no matter the energy value of the new configuration. In contrast with the results of (a) and 
(b), in this case we find that the parabolas become steeper as F — > F^, corresponding to an 
increase in A as the depinning transition is approached as in Eq. (^) p!8|^0|. 

A typical system to which this study may be relevant is the motion of a single fiux line 
in a type-II superconductor, directed along the external magnetic field H . At moderate 
fields, when the separation of the vortex lines is sufficiently large, the intervortex interaction 
can be neglected. In this regime, the dynamics of the vortex phase can be understood by 
studying the motion of a single vortex. We argue that the condition < i^y can be met in 
some anisotropic superconductors. Thus our results might be important in understanding 
the driven diffusion of the flux line. The variation of the velocity with tilt suggests that a 
tilted external magnetic fleld H changes the velocity of the flux line, the effect decreasing 
as we approach the depinning transition. 

In summary, we present a model to describe the motion of an elastic anisotropic string 
in a disordered medium. We flnd that when transverse fluctuations (along the driving force) 
are energetically more favorable than longitudinal fluctuations, the string is described by 
Eq. (|^) not only at the depinning transition but only in the moving phase. However, if 
longitudinal fluctuations are more favorable, a kinematic nonlinear term is induced so that 
its coefficient vanishes as we approach the depinning transition. The directed percolation 
depinning universality class is obtained when a RSOS condition is applied that favors the 
growth of regions with large local slopes. This last result is shown to be vahd for any value 
of the anisotropic parameter e. 

We thank L. A. N. Amaral, R. Cuerno, K. L. Lauritsen, and S. Tomassone for valuable 
discussions. The Center for Polymer Studies is supported by the NSF. 
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FIGURES 

FIG. 1. The discrete version of the elastic string is composed of L massless beads interacting 

via springs. A driving external force F acts in the y direction. Point-like quenched disorder (not 
shown) is introduced at each site on the lattice. The beads are allowed to move in the x and y 
direction and therefore they develop overhangs. 

FIG. 2. Plot of average velocity versus the average tilt of the string for different values of the 
reduced force ranging from / = 0.03 (bottom curves) to / = 0.20 (top curves). Results are averaged 
over 200 independent realizations of the disorder. The system size is L = 250 and the strength of 
the disorder is (5 = 3. (a) The elastic constants are i/^ = 0.1 and Vy = 1 {e = 10). The opening of 
the parabolas as the depinning transition is approached indicates that a nonlinear term is present 
in the equation of motion, and that its value converges to zero at the depinning transition. The 
continuous lines are the best polynomial fits to the curves. We note that the observed parabolic 
dependence cannot be a lattice effect since in this case one expect f(m = 0) = v{m = 1). (b) 
The same plot for the case when longitudinal fluctuations are not energetically favorable, so that 
overhangs are not observed. The elastic constants are = 1.0 and Uy = 0.1 (e = 0.1). The 
horizontal lines indicate that the velocity is independent of the average tilt of the string, (c) Tilt 
dependence of the average velocity of the string when the RSOS condition is applied showing the 
closing of the parabolas indicating a diverging A term at the depinning transition. Here is shown 
the isotropic case = Vy = 1 (e = 1), although the divergency of A is shown to be independent of 
the parameter e. 
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